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The spin relaxation time T\ in zinc blende GaN quantum dot is investigated for different magnetic 
field, well width and quantum dot diameter. The spin relaxation caused by the two most important 
spin relaxation mechanisms in zinc blende semiconductor quantum dots, i.e. the electron-phonon 
scattering in conjunction with the Dresselhaus spin-orbit coupling and the second-order process of 
the hyperfine interaction combined with the electron-phonon scattering, are systematically studied. 
The relative importance of the two mechanisms are compared in detail under different conditions. It 
is found that due to the small spin orbit coupling in GaN, the spin relaxation caused by the second- 
order process of the hyperfine interaction combined with the electron-phonon scattering plays much 
more important role than it does in the quantum dot with narrower band gap and larger spin-orbit 
coupling, such as GaAs and InAs. 

PACS numbers: 73.21. La, 71.70.Ej, 85.75.-d 



I. INTRODUCTION 

The wide-band-gap group III nitride semiconductor 
GaN has emerged as a leading material for a variety of 
new devices, ^'^ ranging from the blue laser~i4 to high- 
power electronic devices,^ by utilizing its electronic and 
optical properties. Recently the magnetic properties of 
GaN-based nanostructures have also attracted much at- 
tention, due to the potential application in spintronic 
device.^ Understanding the carrier spin relaxation mech- 
anism in GaN is of great importance in the design and 
the realization of GaN-based spin device. So far, much 
effort has been devoted to the experimental study of 
the spin relaxation in different GaN structures, includ- 
ing GaN epilayers^iSii^'^" GaN quantum wellsiiii^ii^ and 
GaN quantum dots (QDs).^"' Most of these works fo- 
cus on the spin life time in the hexagonal wurtzite 
GaN structures, which are easier to grow than the cu- 
bic structures. However, the spin-orbit coupling (SOC) 
in wurtzite GaN structure is much larger than that of 
cubic GaN due to the strong built-in electric field caused 
by the spontaneous and piezoelectric polarizations . 
The electron/exciton spin life time of different wurtzite 
GaN nanostructures ranges from a few to a few hundred 
picoseconds. While the exciton spin relaxation time is of 
nanoseconds for in the cubic GaN epilayer^^ and is even 
longer in cubic GaN QD.^"* On the theoretical side, spin 
relaxation times of electron and hole in bulk cubic GaN 
are calculated and are found to be two or three orders 
of magnitude longer than those in GaAsJ^ii^ However, 
the electron spin properties in cubic GaN QDs are less 
well understood and many questions, such as what the 
dominant spin relaxation mechanism is, remain open. In 
this paper, we will systematically study the electron spin 
relaxation in cubic GaN QD under different conditions. 

There are many spin relaxation mechanisms in 
QDs;iSi2£i^i2^i2^i2i In cubic semiconductor QDs, the 
most important two mechanisms are: (1) the electron- 



phonon scattering in conjunction with the SOC, and 
(2) the second-order process of the hyperfine interaction 
combined with the electron-phonon scatteringj^Sii^^iSi In 
GaAs QD, it was shown that the first mechanism is the 
dominant spin relaxation mechanism for quite wide range 
of parameters due to the large SOC.^- Since the SOC in 
GaN is much smaller than that of GaAs, which of these 
two mechanisms dominates spin relaxation need to be 
further examined. 

We organize the paper as following. In Sec. [H] we set 
up the model and give the Hamiltonian. The two most 
important electron spin relaxation mechanisms are dis- 
cussed and the formula of the corresponding spin relax- 
ation rates are presented. We then calculate the spin 
relaxation rates of a QD embedded in a narrow quantum 
well analytically using perturbation theory in Sec. IIIII 
We further present the exact spin relaxation rates under 
different conditions by numerical method in Sec. llVl and 
summarize in Sec. [V] 

II. MODEL AND SPIN RELAXATION RATE 

We consider one electron spin in a single GaN QD em- 
bedded in a quantum well with well width a. A magnetic 
field B is applied. The Hamiltonian of the system com- 
posed of the electron and the lattice is given by: 

Ht = H, + Hl + H,l , (1) 

where He , Hl and H^l are the Hamiltonians of the elec- 
tron, the lattice and their interaction respectively. The 
electron Hamiltonian He can be written as: 

He = Ho + Hso 

- [^ + Veix,y) + V,iz) + Hz]+Hso , (2) 
2m* 

where Hq is the Hamiltonian without the SOC, m* is 
the electron effective mass and P = —ih\/ -I- - A is the 
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kinetic momentum with A = B x r. Vz{z) is the quan- 
tum weh confinement. In this paper, it is assumed to 
be a hard wall confinement with width a. Vc{x,y) ~ 
^m*ujQ{x'^ + y^) is the in-plane confinement of QD with 

diameter do — yhnfrnFuM^. Hz = ^g/isB • a is the 
Zeeman energy with [Ib and a being the factor of 
electron, Bohr magneton and Pauli matrix, respectively. 
Hso is the Hamiltonian of the SOC. In cubic GaN the 
dominant SOC term is Dresselhaus term,-2S. which reads^^ 
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Uso = J^loiiPyPxPy - PzPxPz)a^ + {PzPyPz 
PxPyPx)^y ^" {^PxPzPx ~^ PyPzPy)^ z\ 



(3) 



with 7o being the Dresselhaus coefficient. The eigen 
wave function \t) and the eigen energy = 1, 2, • • • ) 
of Hf. can be obtained from the perturbation theory 
or from the exact diagonalization method, using the 
eigenstates of Hq as basis. The Hamiltonian of the 
lattice is consisted of two parts: H]^ = Hp^ + Hnuciei- 
Hph = X]qi) ^qifflqrjCtqi; rcprcseuts the Hamiltonian 
of the phonons with u^r) standing for the phonon 
energy spectrum of branch rj and momentum q 
and <iLJaqr) being the corresponding phonon cre- 



ation/annihilation operator. Hn 



E, 7/B-I, is 



the Zeeman term of the the lattice nuclear spins in 
the external magnetic field with 7/ and Ij denoting 
the gyro-magnetic ratio and spin of j-th nucleus re- 
spectively. The interaction between the electron and 



the lattice also has two parts: HeL — H^p + Hei- 
Hep is the electron-phonon scattering and is given by 

Hep = Eq„^V)( 



qj; ^"qifV'^qr) 



-qr; 



)e^'^'^, where Mq^ is the 
matrix element of the electron-phonon interaction. 
l-^qsiP = KE?q/2pVsi for the electron-phonon coupling 
due to the deformation potential. For the piezoelectric 
coupling, |Mqp/p = {■i2hTT'^e'^el^/ pVsi)[{iqxqyqzY / q'] 
for the longitudinal phonon mode and X]i=i 2 l-^qptj P ~ 
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qtqi 
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{iqxqyqz)^ / q^] for the two transverse modes. Here 
E! stands for the acoustic deformation potential; p is the 
GaAs volume density; 614 is the piezoelectric constant 
and K denotes the static dielectric constant. The 
acoustic phonon spectra Wq/ = Vgiq for the longitudinal 
mode and Wqt = Vstq for the transverse modes with Vsi 
and Vst representing the corresponding sound velocities. 
Hei is the electron- nucleus hyperfine interaction H^i, 
which can be written as H^i — ^"^oS ■ Ij<^(r — Rj); 
where wq is the volume of the unit cell of the lattice, S 
is the spin of the electron, r and Rj are the position of 
the electron and the j-th nucleus, respectively. A stands 
for the hyperfine interaction coupling constant. 

In the Hamiltonian [Eq. ([1])], we only include the terms 
related to the two dominant spin relaxation mechanisms. 
One is the electron-phonon scattering in conjunction with 
the Dresselhaus SOC. The SOC mixes the spin-up and 
-down states to form the majority spin-up and -down 
states. The direct coupling to the phonon causes the 
transition between the majority spin-up and -down states 



and results in the spin relaxation. The transfer matrix el- 
ement is Mqri ■ This spin mechanism is referred as "Mech- 
anism I" hereafter. The other is the second-order process 
of the hyperfine interaction combined with the electron- 
phonon interaction in which not only the SOC mixes the 
spin- up and spin-down states, but also the nuclei flip the 
electron spin. As the phonon compensates the energy 
difference, this mechanism also leads to spin relaxation. 
In the following, it is called "Mechanism H" . The transfer 
matrix between states |£i) and [^2) of Mechanism II can 
be written as 



el—ph 



E 



{£2\Hep\m){m\Hei\£i) 



E 



£li — Sm 

{e2\Hei\m){m\Hep\ei) 



££2 ~ ^r, 



(^1 



(4) 



with 



Xq,, = \£2 



E 



(^2|Mq^e^i"-|m)(m|i7e/|^i) 



■E 



{e2\H,i\m){m\M, 



qjji: 



,iq-r| 



(^il ,(5) 



where the summation of |m) runs over all possible inter- 
mediate states. 

To calculate the spin relaxation time, one can use 
the perturbative approach based on the calculation of 
the transition rates from Fermi's golden ruleJ^iSSii^'^iiSl 
Non-perturbative calculation using equation of motion 
method has also been proposed to study the spin relax- 
ation of the system with large SOC at high temperature 
regime.^^'^** For the system with weak SOC at low tem- 
perature regime, these two approaches produce the same 
results. In the cubic GaN QD, since the SOC is pretty 
weak)^ the perturbative approach gives sufficient accu- 
rate spin relaxation rate and is therefore adopted in the 
present work. 

Using the Fermi's golden rule, one can obtain the spin 
relaxation rate asi^ 

Tr' = E(/^+r»+-/- + . (6) 

Here fi± is the Maxwell distribution since we study the 
spin relaxation of single electron confined in the QD. 
"+/—" stand for the states with the majority up/down- 
spin. The scattering rate Ti^f reads: 

27]- 

Ei^/ = — ^|(/|A'q^|i)|^[nq^(5(e/-e, -^iwq^) 



qr) 



+ l)S{ef ~ Si + hujq,,)] , 



(7) 



where Xq,^ = Mq^^e'^'^'^ and X^,^ = Mqij for Mechanism I 
and II, respectively, nq^ is the Bose distribution function 
for phonons. 
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III. ANALYTICAL RESULTS 

Before presenting the full exact diagonalization result, 
let us first look at the analytical result of the spin re- 
laxation rate of a QD embedded in a narrow quantum 
well by perturbatively solving the electron Hamiltonian 
[Eq. ^] to the second order of the SOC. 

Due to the symmetry of the QD in the x-y plane, B 
can be assumed to be {B^^,0, B±) with _B|| = Bsm9 
and B± = Bcos9 being the components along the a;- 
and z-axis, and 9 representing the angle between the 
magnetic field direction and the z-axis. The eigen- 
state of the electron Hamiltonian without the SOC {Hq) 



luzula) is characterized by the quantum number of quan- 
tum well confinement, radial, angular and spin freedoms 
1,2,---), n (= 0,1,---), I 0,±1,---) and 

a (= ±1) respectively, whose energy is -B„^„zo- = Im'a'^ + 

{2n+ \l\ + l)hn + lhujBs_ +(tEb, with fl = ^Juj^ + , 

u!b — eB±/{2m*) and Eb — ^qubB. In the narrow 
quantum well (do ^ a), the distance of different states 
is so large that only the lowest state is relevant. Un- 
der this approximation, the spin orbit coupling can be 
expressed as Hgo = j-Jo{'^/0')'^{—PxO'x + PyO'y)- Up to 
the first order perturbation, the lowest two eigenstates 
of the electron with the SOC are 



l^i) = 1001) -HS^ltooi 1011) -^?ol^,ooil01-l)+Stn,ooilO- 11) +^?oL"^^^^^ , 

= |00-l)+6o'itoo-il011)-^^ol-i,oo-il01-l)-So^n,oo-ilO-ll)-S^^^^^^ (8) 

where 6;^'±„,,,^, = qhjca{'^±ujB^ /il) sin 9 /{Enla-En>l, a'), 13li%,^,i,^, = i^hjM^±^B^/^){'^+cos9)/{Enia~En'i'a') 
and ^S^jct n';v' ^ i^fi'ycOi{l±ojB± /^){^ — cos9) /{Enia — En' I' a')- The corresponding eigen energies of these states read 

£1 = -Eooi + |^ditooiP(^oii - -Booi) - |Soi-i.ooiP(-^oi-i - ^^ooi) + l^d-ii,ooi P(-^o-ii - £^001) 

+ 1^0-1-1,001 P(-^o-i-i --B001) , 
£2 = £^00-1 + l^on,oo-iP(-^oii - £00-1) - |^di-i.oo-iP(-^oi-i - £00-1) - I^o-ii,oq-iP(^o-ii - -^oo-i) 

-|^(lli_i,oo-iP(^o-i-i - £00-1) ■ (9) 

I 



It is noted that in the above equations, we have included 
the second-order correction of the SOC on the energy, 
which is crucial to the study of the spin relaxation using 
perturbation method as pointed out by Cheng et al?^ 
It is also noted that the state index is dropped in 
the above equations since it is always 1 in the narrow 
quantum well approximation. and [^2) are the low- 
est majority spin-up and spin-down states respectively. 



At the low temperature regime, the electron mainly dis- 
tributes on these two states. Therefore, Ti basically 
equals the spin relaxation time between these two states. 
It is further noted that at low temperature regime, the 
main electron-phonon scattering comes from the electron 
coupling to the transverse phonon via piezoelectric field. 
With these approximations, the spin relaxation rate due 
to Mechanism I is given by 



Fi = c{2nq + l)aq d9' sin^ 9' {sin'^ 9' + 8 cos'^ 9')e-'^' ^'"^ ^'^^I^{-qaa cos 9')[2P^ 
Jo 2 

+ + Pi - 2Vl)\q^ sin^ 9' + {Vl + Vi + 2Pf sin^ 0'] , (10) 

I 



where c = TTe'^e'l^/{hDv%K'^), q = AE/{hvsta) with 
AE = |£2 - eil and a = ^JmF^lJh. I(x) = 
TT^ sin(a:;)/[x(7r — a;)(7r denotes the form factor along 
z-direction due to the quantum well confinement. In the 
above equation, V\ = A\ + A2 — Az — A4 with Ai — 
1^011, ooi^oli,oo-il' — I'Sgi.! 001^01-1,00-1 1' -^3 — 



^0-11,001^0-11.00-1! ^n-^i — 1^0-1-1,001^0-1-1,00-1 1; 

"^2 = ^I'^oii.oo-il + 1^0-11, oo-il ^ I'^oi-i.ooil + 

1^0-1-1,001 1' "^3 ^ ^l^oll,00-ll ^ 1^0-11, OO-ll + 
1^01-1,001 1 + 1^0-1-1,001 1' ^4 = Ci - C2 + C3 - C4, 

7^5 = Ci - C2 - C3 -I- C4 with Ci = |Sdi7^ooi^o-^ii,oo-iL 
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^2 — 1^01-1, 001'^0-l-l,00-ll' ^3 — 1^0-11,001^011,00-1 1 

and C4 = \Bq'_i_i ooi'^oi-i 00-1 1- Using the material pa- 
rameters of GaN QD and in consideration of the relative 
small magnetic field, one can write down the spin relax- 
ation rate due to this mechanism at zero temperature for 
relative small dot: 

Ti cx a'-^dlB^ (l + cos^ 0) , (11) 

which indicates that for fixed magnetic field magnitude, 
the spin relaxation under the perpendicular magnetic 
field is two times of that under the parallel magnetic filed. 
It should be noted that for 9 = case (the magnetic field 
is along the z-axis), to the leading term, the magnetic 
field dependence of Fi obtained here is in accordance 
with that obtained in Refs. [l^ and[20l. By assuming that 
the nuclei spins are independent to each other and are in 
equilibrium state, the spin relaxation between \£i) and 
{£2) induced by Mechanism II, with the mediation of the 
lowest available state, can be written as 

T2 - ( ^ fliI+l)voa^a-^c{2nq + l)q^ 

£2 — £3 

null 

X / d6l'sin^6l'(sin4 6l'-^8cos''6l') (12) 

xe-«'^'"''''/2/2(igaacos0') , 

which at zero temperature gives 

Ta oc C3a-^4B^ . (13) 

The ratio of the spin relaxations due to these two mech- 
anisms is therefore 

Ti/Fa cx a-^4B^ , (14) 

which gives a guideline to determine which mechanism is 
more important at different conditions. It is therefore ex- 
pected that Mechanism II is more important for smaller 
QD embed in wider quantum well under weaker magnetic 
field. 



IV. NUMERICAL RESULTS 



TABLE I: Parameters used in the calculation. 



p 


6.095 X 10^ kg/m^ 




8.5 


Vst 


2.68 X 10^ m/s 


9 


2.06 


Vsl 


6.56 X 10^ m/s 




8.3 eV 


ei4 


4.3 X 10" V/m 


m* 


0.15mo 


A 


45 fieV 


I 
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The perturbation method gives qualitative results for 
us to understand the overall behavior of spin relaxation in 
GaN under different conditions. However, in the pertur- 
bation calculation, states with higher energy are dropped 



to get a manageable analytical result. It should be noted 
that, for the spin relaxation caused by Mechanism II, 
the contributions of higher intermediate states and the 
lowest one are of the same order in regard to hyper- 
fine interaction strength. Moreover, for QD embedded 
in wider quantum wells, contribution of the higher Uz 
states to the spin orbit coupling can not be neglected. 
It is expected that for c?o ^ the spin relaxation due 
to Mechanism I can be different from the perturbative 
results. It is therefore necessary to check the accuracy of 
the perturbative calculation by comparing to the exact 
diagonalization with sufficient basis functions included. 

In Fig. [1] we present the spin relaxation rates as func- 
tions of well width in GaN QD under different condi- 
tions obtained by the exact diagonalization and pertur- 
bation. The material parameters of GaN are listed in Ta- 
ble HJ^i^i^i^ The Dresselhaus coefficient 70 is chosen to 
be 0.51 A3-eV according to the latest calculation.— It is 
seen that the perturbation results describe the qualitative 
behavior of the spin relaxation pretty well. For the cases 
we study here, the spin relaxation caused by the electron- 
phonon scattering in conjunction with the SOC from the 
perturbation is very close to the exact diagonalization 
result in narrow quantum well. When the well width 
becomes larger, the difference between perturbative and 
exact diagonalization result also grows as contribution of 
the higher states becomes more and more important. 
On the other hand, for the spin relaxation caused by the 
hyperfine interaction together with electron-phonon in- 
teraction, the difference between perturbative and exact 
diagonalization results almost doest not change with the 
well width. For this mechanism, the relaxation rate from 
exact diagonalization method is always about one order 
of magnitude higher than that obtained from the pertur- 
bation calculation for the cases we study. This indicates 
the contribution of the higher states are important to 
the quantitative calculation of the spin relaxation. In 
the following, we only present the results of exact diag- 
onalization unless otherwise specified. We now focus on 
how the spin relaxation induced by the two mechanisms 
change with a. It is seen that the spin relaxations in- 
duced by the two mechanisms both decrease with a. The 
spin relaxation due to Mechanism I decreases much faster 
than that due to Mechanism II. As a result, Mechanism 
II becomes more and more important when the quantum 
well width increases. This can be understood from the 
perturbation result. As one can see from Eqs. (jlip and 

that the relaxation rate Fi decreases with a as 
while F2 is proportional to a~^. We then pay attention 
to the relative importance of these two mechanisms. For 
the vertical-magnetic- field case in Fig.[TJa), for S = 5 T, 
spin relaxation due to Mechanism I is always the dom- 
inant spin relaxation mechanism. When B decreases to 
0.5 T, Mechanism II almost dominates the spin relax- 
ation except at very small well width (a < 4 nm). For 
the parallel-magnetic-field case in Fig. [lib). Mechanism 
II is even more important and dominates the spin relax- 
ation for a > 6 nm and > 2.5 nm when _B = 5 T and 0.5 
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T, respectively. This is quite different from the cubic ma- 
terials with narrower band gap and larger SOC such as 
GaAs, in which the spin relaxation due to Mechanism I is 
usually 2-3 orders of magnitude stronger than that due to 
Mechanism II. But thanks to the small SOC, the spin re- 
laxation caused by the nuclei plays much more important 
role in GaN QD. It is also worth noting that the hyperfine 
interaction and the SOC can also cause spin dephasing. 
Previous studies on GaAs QD have shown that the hy- 
perfine interaction usually dominates the spin dephasing 
at low temperaturci^i^ It is expected that the spin de- 
phasing in GaN QD is also dominated by the hyperfine 
interaction due to the very small SOC in this material. 
Our numerical results using the approach in Ref . show 
that this is indeed true, e.g., for QD of a = 5 nm and 
do = 10 nm, T2 induced by the hyperfine interaction is 
about five orders of magnitude shorter than that induced 
by the SOC under parallel magnetic field of 0.5 T when 
T = K. As we are interested in the difference between 
GaN and GaAs QDs, we will not further address the spin 
dephasing in the paper. 
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FIG. 1: (Color online) Spin relaxation rate vs. the well width 
a in the presence of (a) perpendicular and (b) parallel mag- 
netic fields with B = 0.5 T (solid curves) and _B = 5 T (dot- 
ted curves). In the calculation, do = 10 nm. Black (dark) 
curves — exact diagonalization results; Red (light) curves — 
perturbation results. Curves with v — Tj~^ induced by the 
electron-phonon scattering in conjunction with the SOC; 
Curves with • — T^^ induced by the second-order process of 
the hyperfine interaction combined with the electron-phonon 
scattering. 

In Fig. [5] the QD diameter dependence of the spin re- 
laxation is presented under the magnetic field perpen- 
dicular (a) and parallel (b) to the well plane. Both re- 
laxation rates increase with the increase of dot size but 
with different speeds: Fi oc dp and F2 oc d^. As a re- 
sult. Mechanism I becomes more important as the size 
of QD grows. One can see from Fig. [2] that, under the 
low magnetic field {B = 0.5 T) we show here, Mecha- 
nism II plays a very important role, or even dominates 
the spin relaxation for all QD whose diameter is smaller 
than 11 nm. 

In Fig. [Sja) and (b) the spin relaxation rates induced 
by the two mechanisms are plotted as functions of the 
perpendicular and parallel magnetic fields respectively. 
In each figure, the results are shown for both narrow well 
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FIG. 2: Spin relaxation v&. the QD diameter do in the pres- 
ence of (a) perpendicular and (b) parallel magnetic fields at 
two well widths: a = 5 nm (solid curves) and a = 10 nm 
(dotted curves). In the calculation, B = 0.5 T. Curves with 
V — Tj~^ induced by the electron-phonon scattering in con- 
junction with the SOC; Curves with • — T^^ induced by the 
second-order process of the hyperfine interaction combined 
with the electron-phonon scattering. 



[a — h nm) and relatively wide well (a — 10 nm). It is 
noticed that the effect of each mechanism increases with 
the magnetic field as predicated by Eqs. (fTTj) and (fT^ . 
Then we pay attention to the relative importance of the 
two mechanisms. When the magnetic field is along the 
2;-direction, it is seen from Fig. [3] that Mechanism I is 
dominant when large vertical magnetic (B > 0.5 T) is 
applied. However, when the magnetic is along x-axis, for 
small well width (a = 5 nm), Mechanism I is dominant 
for large magnetic field. For wider quantum well (a = 10 
nm), Mechanism II dominates the spin relaxation when 
B < 2.5 T and is comparable to Mechanism I for larger 
magnetic field. 

We then turn to study how the direction of the applied 
magnetic field changes the spin relaxation. In Fig. 31 
we show the spin relaxation rates as functions of the 
angle Q between B and the z-axis for a fixed magnetic 
field amplitude. It is seen that these two mechanisms 
depend on the direction of the magnetic field quite dif- 
ferently. The spin relaxation induced by Mechanism I 
has a maximum when the magnetic field is along the z- 
direction. With the increase of 6*, it decreases gradually 
and reaches the minimum when B is in the x-y plane. On 
the other hand, the spin relaxation induced by Mecha- 
nism II almost keeps unchanged with Q. This can be 
understood from the perturbation result. As we can see 
from Eqs. (|lip and (fT3)) that the relaxation rate Fi con- 
tains the term of (1 -I- cos^ 6*), which has the largest value 
for 6* = and the smallest value for Q = 7r/2 for the 
condition we considered, while r2 is almost independent 
of d. Overall, the changes of the spin relaxation rates 
in GaN QD are mild when the direction of the magnetic 
field changes for both mechanisms. This is quite differ- 
ent from that in GaAs QD, where the spin relaxation in- 
duced by Mechanism I with the perpendicular magnetic 
field can be several orders of magnitude larger than that 
with the parallel magnetic field, This is because in 
GaAs material, the SOC is usually comparable or even 
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larger than Zeeman splitting and therefore the magnetic 
field direction changes the eigen energy remarkably. Con- 
sequently, the difi^erence between the maximum and the 
minimum of the spin relaxation rates induced by Mecha- 
nism I can be several orders of magnitude different when 
the direction of the magnetic field changes. However, due 
to the small SOC in GaN, the energy difference between 
the lowest two eigen states is determined by the Zeeman 
splitting and therefore the change of the spin relaxation 
with the magnetic field direction is much milder. 
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FIG. 3: Spin relaxation vs. (a) perpendicular and (b) parallel 
magnetic field at two well widths: a = 5 nm (solid curves) 
and a = 10 nm (dotted curves). In the calculation, do = 10 
nm. Curves with v — induced by the electron-phonon 
scattering in conjunction with the SOC; Curves with • — T^^ 
induced by the second-order process of the hyperfine interac- 
tion combined with the electron-phonon scattering. 
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FIG. 4: Spin relaxations vs. 6. In the calculation, a = 5 nm, 
B — 0.5 T and do = 10 nm. Curve with v — induced by 
the electron-phonon scattering in conjunction with the SOC; 
Curve with • — T{'^ induced by the second-order process of 
the hyperfine interaction combined with the electron-phonon 
scattering. Note the scale of the spin relaxation induced by 
the second mechanism is at the right hand side of the frame. 

We further investigate how the spin relaxation changes 
with the temperature. The results are shown in Fig. [5] 
One can see that spin relaxations induced by the two 
mechanisms both increase with the temperature. For 
low temperature regime, the relative importance of each 
mechanism remains unchanged. That is. Mechanism I is 
more important when the magnetic field is perpendicu- 
lar to the well, while Mechanism II usually plays more 
important role for the parallel magnetic field. Both are 
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FIG. 5: (Color online) Spin relaxation vs. temperature T 
in the presence of (a) perpendicular and (b) parallel mag- 
netic fields. In the calculation, a = 5 nm, do = 10 nm and 
B = 0.5 T. Black (dark) curves — exact diagonalization re- 
sults; Red (light) curves — perturbation results. Curves with 
V — induced by the electron-phonon scattering in con- 
junction with the SOC; curves with • — T^^ induced by the 
second-order process of the hyperfine interaction combined 
with the electron-phonon scattering. 

approximately proportional to [2nq(r) -f 1] which is con- 
sistent with the perturbative results, i.e., Eqs. pUj) and 
(I13p . However, when the temperature rises high enough 
(T > 10 K), the spin relaxation induced by Mechanism 
I increases much quicker than Mechanism II. For both 
parallel and perpendicular magnetic fields, Mechanism II 
dominates the spin relaxation for low temperature while 
Mechanism I has larger contribution for high tempera- 
ture. In order to understand the different temperature 
dependences of relaxations, we also show the spin re- 
laxation rates from perturbation calculation in the same 
figure. It is seen that the perturbation result and the 
exact diagonalization result of the spin relaxation due to 
Mechanism II agree with each other qualitatively in the 
temperature regime we study. However, the spin relax- 
ation of exact diagonalization due to Mechanism I de- 
parts from the perturbation prediction in the high tem- 
perature regime. This indicates that the perturbation 
method is no longer reliable for Mechanism I and should 
not be used to obtain the spin relaxation rate under high 
temperature. This is understandable, because at low 
temperature regime, the electron distribution at the high 
levels is negligible and only the lowest two Zeeman split- 
ting levels are involved due to the small SOC. Therefore 
perturbative method is accurate enough to study the spin 
relaxation caused by Mechanism I. With the increase of 
the temperature, electron can occupy higher energy levels 
with larger SOC. As a result, the perturbation method is 
no longer adequate to study the transition rates induced 
by Mechanism I. 



V. CONCLUSION 

In conclusion, we have investigated the spin relax- 
ation time Ti in cubic GaN QD under different condi- 
tions by the perturbation and exact diagonalization ap- 
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proaches. Two leading spin relaxation mechanisms, i.e., 
the electron-phonon scattering in conjunction with the 
SOC and the second-order process of the hyperfine inter- 
action combined with the electron-phonon scattering, are 
considered. We systematically study how the spin relax- 
ations induced by the two mechanisms change with the 
well width a, magnetic field B and quantum dot diame- 
ter do- Our results show that, the ratio of these two spin 
relaxation rates is proportional to a^^B^dg in the low 
temperature regime when the quantum well constraint 
is strong enough. Due to the small SOC, the spin relax- 
ation caused by the second-order process of the hyperfine 
interaction combined with the electron-phonon scatter- 
ing plays much more important role in GaN material. 
Only when the well width a is small enough and/or the 
magnetic field B and QD diameter do are large enough, 
the electron-phonon scattering in conjunction with the 
SOC may dominate. Furthermore, how the direction of 
the applied magnetic field changes the spin relaxation 
are investigated. The spin relaxation induced by the 
electron-phonon scattering in conjunction with the SOC 
has a maximum when the magnetic field is along the z- 



direction and reaches the minimum when the magnetic 
field is in the x-y plane. Nevertheless, the spin relaxation 
induced by the second-order process of the hyperfine in- 
teraction combined with the electron-phonon scattering 
keeps almost unchanged with the magnetic field direc- 
tion. We also discuss the temperature dependence of the 
spin relaxation due to the two mechanisms. At high tem- 
peratures, the spin relaxation induced by the electron- 
phonon scattering in conjunction with the SOC is always 
dominant. 
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